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Abstract 

We develop a formalism to carry out coarse-grainings in quantum field theo- 
retical systems by using a time-dependent projection operator in the Heisen- 
berg picture. A systematic perturbative expansion with respect to the in- 
teraction part of the Hamiltonian is given, and a Langevin-type equation 
without a time-convolution integral term is obtained. This method is applied 
to a quantum field theoretical model, and coupled transport equations are 
derived. 

I. INTRODUCTION 

Obtaining descriptions of nonequilibrium processes in quantum field theory is an impor- 
tant problem and has been actively studied. In the treatment of such systems, it is often 
important to carry out reductions or coarse-grainings of irrelevant degrees of freedom. In 
classical systems, it is widely believed that a system of macroscopic size composed of many 
microscopic variables exhibits rather simple macroscopic behavior that can be described in 
terms of only a few macroscopic variables. It is therefore not unreasonable to conjecture that 
such an elimination of the irrelevant information can be developed for application to quan- 
tum systems and help in the analysis of nonequilibrium quantum processes. The projection 
operator method is one well-known method for carrying out coarse-grainings systematically. 
|T|-[5]] In this method, after the elimination of irrelevant information by means of a projection 
operator, a kind of master equation in the Schrodinger picture or a Langevin-type equation 
in the Heisenberg picture is obtained. Famous examples are the Nakajima-Zwanzig and 
Mori equations. Recently, the unification and generalization of the treatments used in 
these two picture has been realized. 

In the derivation of kinetic equations from the microscopic point of view, we usually 
assume that the relevant part of the system interacts with an irrelevant subsystem, which is 
regarded as a heat reservoir in thermal equilibrium. || However, in a general nonequilibrium 
system, we can expect the irrelevant subsystem to have some time dependence. Further- 
more, it is often convenient to have a description in which two or more coupled systems are 
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considered to be on an equal footing. In particular, in quantum field theory, both relevant 
and irrelevant parts of the system may be composed of infinite numbers of degrees of free- 
dom, and therefore, it is not clear whether we can regard one subsystem as a reservoir for 
the other or not. To elucidate this situation, it is convenient to introduce a time- dependent 
projection operator. 

The time-dependent projection operator was first applied by Robertson. He attempted 
to project the complex behavior of the relevant density matrix onto a local equilibrium den- 
sity matrix and derived a master equation. Ochiai attempted to obtain a kinetic equation by 
using a similar projection. || Robertson's time-dependent projection operator was improved 
by Kawasaki et al. ffl Willis et al. considered coupled systems and introduced another time- 



dependent projection operator. [Tjj A similar projection was implemented by Grabert et 



al. jn] Shibata et al. derived a systematic perturbative expansion formula for a master 
equation without time-convolution integral terms. [12[ 



In all of the above cited works, formulations were developed in the Schrodinger picture, 
and master equations were derived. Langevin-type equations are obtained when we apply 
the projection operator method in the Heisenberg picture. Such equations were studied by 



Grabert (13| and Furukawa. |L4| They carried out an exact transformation of the Heisenberg 
equation of motion to derive a Langevin-type equation. However, to apply this approach 
to concrete phenomena, such an exact formulation is not convenient, and we must resort 
to perturbative calculations. In this paper, we develop a systematic perturbative expansion 
formula of a Langevin-type equation. It is known that both equations with and without a 
time-convolution integral term can be derived using the projection operator method. @-[5| 
In this study, we investigate an equation without a time-convolution integral term. Further- 
more, to investigate the validity of our formalism, we apply it to a quantum field theoretical 
model and attempt to derive coupled transport equations. 

This paper is organized as follows. In §HJ, the projection operator method with a time- 
dependent projection operator is developed. In §|TJ, we apply our formalism to a quantum 
field theoretical model that is composed of two bosons. We assume local equilibrium and 
derive coupled transport equations. A summary and conclusions are given in § |PVl . 

II. TIME-DEPENDENT PROJECTION OPERATOR METHOD 

Our starting point is the Heisenberg equation of motion, 

±0{t)=i[H,0®) (!) 
= iLO{t) (2) 
— > 0(t) = e iL{t - t0) O(t ) } (3) 

where L is the Liouville operator and t Q is the time at which we prepare an initial state. The 
Heisenberg equation contains complete information of the time-evolution of the operator, 
but in general, it is difficult to solve exactly when there are interactions. For this reason, it 
is necessary to carry out a reduction or a coarse-graining of the irrelevant information. For 
this purpose, we introduce the time-dependent projection operators P(t) and Q(t), which 
are related as 
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Q{t) = l-P(t). (4) 

The projection operator P(t) helps us to project any operator onto the P-space, which con- 
sists of the relevant degrees of freedom. From Eq. (|3|), one can see that the time dependence 
of the operators is determined by e tL ^~ to \ This yields 

dt 

= e lL{t - t0 \P(t) + Q(t))tL. (5) 
From this equation, we can derive the two equations 

d lL (t-t ) p ^ = e iL(t-t )( p ^ + Q(t))iLP(t) + e lL{t - to) P(t), (6) 

(JjL 

diL(t-t )Q^ = e ii(t-t )(p( t ) + Q( t ))iLQ(t) + e lL{t ' to) Q(t), (7) 

(aiL 

where P(t) = dP{t)/dt and Q{t) = dQ{t)/dt. Equation (g) can be solved for e iL{t ~ to) Q: 
e iL ^Q(t) = Q(t )S dSLQiS) + fdse^HQ{s)+P{sYLQ{s))eS dTLQ(T) 



i I dsLQ(s) f . 

Q(t )e> + e lL ^- t0 \P(t) + Q(t))Z(t,t, 

dsLQ(s) . _ , . I 

Wt 0) e^ +e ^-to)p (OS(t)to)} ___ ) (9) 



(8) 



where 



£(Mo) = f dse- iL ^{Q(s) + P( S )*LQ( S )}etf drLQ(r) . (10) 

•/to 

i I dsLQ(s) 

Here, the time ordered operator e-/° is defined as 

i f* dsLQ(s) °° /■* /•*! /■*"-! 

J in v ' -i X ^ -nil,! 1 1 I 



1 + Y,i n dh dt 2 --- dt n LQ{t n )LQ{t n ^) ■ ■ ■ LQ{h) 

1 J tQ J Iq J tQ 



n=l 



Note that there is a term including Q(s) that is not observed in the the time-independent 
projection operator method. 

Substituting Eq. (|) into Eq. (|5|) and operating with O(t ) from the right, we obtain 

jO(t) = e iL{t ~ to) P(t)iLO(t ) 

+e iL ^P(t)X(t,t ) l _ * to) iLO(t ) 

i f* dsLQ(s) 1 

+g(*o)e> j—— - iL(t,t )O(t ). (12) 
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This equation has no time-convolution integral. This point is demonstrated at the end of this 
section. We call this the "time-convolutionless" (TCL) equation. When the time-dependence 
of the projection operator is ignored, this equation agrees with Eq. (2-13) of Ref. |J. When 
we substitute Eq. (f|) into Eq. (JD, we obtain an equation with a time-convolution integral 
term, which is called the "time-convolution" (TC) equation. However, we discuss only 
the TCL equation in this paper. 

The TCL equation is exactly equivalent to the Heisenberg equation of motion. Usually, 
we carry out the perturbative expansion of the interaction Hamiltonian and retain only the 
lowest order terms, which are often sufficient to describe the dissipation effect. However, 
the above TCL equation is not convenient to carry out the perturbative expansion. Our 
next task is therefore to rewrite the TCL equation. For this purpose, we restrict the nature 
of the projection operator. For the time-independent projection operator P, the condition 
P 2 = P is satisfied. Contrastingly, in the time- dependent case, the order of operation of the 
projection operators is important. In most previous works, the condition P(ti)P(t 2 ) = P{t\) 
was employed. [[T|-|T2|j fl~4|j However, in this paper, we assume the condition 



P(t 1 )P(t 2 )=P(t 2 ), (13) 

because the projection operator that we use in the next section satisfies this condition. From 
this condition, we can derive several relations: 

Q{ti)P{t 2 ) = 0, (14) 
P(t 1 )Q(t 2 ) = P(t 1 )-P(t 2 ) 1 (15) 
Q(ti)Q(h) = Q(h). (16) 

The total Hamiltonian of the system can be divided into two parts, 

H = H (t, t ) + Hj(t, t ), (17) 

where H (t,t ) and Hi(t,to) are the nonperturbative part and the perturbative part of the 
Hamiltonian, respectively. In a general nonequilibrium process, we can consider the case 
in which the mass of a particle changes with time, and therefore, H (t,to) becomes time 
dependent. The corresponding Liouville operators are defined as 

L (t,t )O=[H (t,t ),O], Lj(t,to)0=[Hj(t,to),0}. (18) 

Now, we assume another condition, 

Q{t 2 )L Q {t l ,t )P{t l )=U. (19) 

With this condition, the form of the nonperturbative Hamiltonian affects that of the projec- 
tion operator, and vice versa. This condition means that the nonperturbative Hamiltonian 
does not drive the operator, once it is projected onto the P-space, into the Q-space, which 
is orthogonal to the P-space. 

With the above properties of the projection operators, the function £(t, t ) can be ex- 
pressed as 
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E(t,to) = t dse- lL{t - s \Q(s) + P(s)iLQ(s)}eS dTLQiT) 

Jt 



= Q(t)- e --"-«Q( ( „ )e <' ,Tl0W 

dsLo(s,to)Q(s) 

= to - ty ^MoMMoJTOJ^lMoje^ . (20) 

Here, the operators C(t, to) and V(t,t ) are expressed as 

C(t,t ) = f/ (t,to)e- 4i{t - to) 

= 1 + J2HT dh dt 2 --- dt n L I (t 1 ,t )L I (t 2 ,t )---L I (t n ,t ), 

n= l J to J to J to 



(21) 



i |, dsLQ(s) 
Jt (TT ( J\~ 1 I 



V(t,t ) =e4 t0 ^ (Ug)- l (t,t ) 

OO „f 



1 + ]T i" / dti / 1 dt 2 • • • / " _1 dtJLfih, t )Lf(t 2 , t )--- L?(t n , t ) } 

-—l J to Jt() J to 

(22) 



where 



t ) = ^LjU^it, t ), (23) 
L?(Mo) = ^(Mo^/Wo 9 )" 1 ^), (24) 

U (t,t ) = e-> , (25) 

f/ Q (Mo) =e> . (26) 

It is noteworthy that the term including the time derivative of the projection operator dis- 
appears. The operator C{t, t ) [V(t, t )] is a time [an anti-time] ordered function of Liouville 
operators. (The details of these expressions are given in Appendix [A].) Then, we have 



(27) 



To derive this expression, we have used mathematical induction. || (The detailed derivation 
appears in Appendix [B].) Substituting the above result into Eq. (|T2|), we obtain 

^0(t) = e iL{t - ta) P{t)iLO{to) 

+e iL { t-t o)p{t)u -i^ to)c(t) t )Q(t) i _ 1 U (t, t )iLO(t ) 

i f* dsLQ(s) I 

+Q(to)e^ to r^Uo) iLO{to) - (28) 

This form of the TCL equation is convenient to carry out the perturbative expansion. 
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When we expand -P£(t, to)/ (1 — £q)) up to first order in the interaction Hi, we have 



e iL{t - to) P{t)iLO{t 



e iL ^P{t)U^\t,t Q )Q{t)U Q {t,t Q )iLO{to) 



+e lL(t - t « ) P(t)U \t, t ) / dsU (s, to^Ms, t )U -\s, t )Q(t)U (t, t Q )iLO{t ) 



_ e iL(t-t ) 



P(t)U -\t,t )Q(t) 



X 



f dsU {s, t )tLj{s, t )U -\s, t )Q(t)U {t, t )tLO{t ) 



i f* dsLQ(s) 1 

+ Q(to)ei° T^oj' 10 ^ 

e ,L(t-to)p (t)f/ -l (t;to) p (t)f/o(t;to)2LO(to) 

+e W-to) P (t)Uo 1 (t, t )P(t) f dsU (s, to^LjU^is, t )Q{t)U (t, t )iLO{t ) 



to 



i \l dsLQ(s) 1 

+Q(t )e> i-V{tM) iL ° {h) ' (29) 



It is easily seen that Eq. ( P9|) does not contain a time-convolution integral. If this did 
contain a time-convolution integral, the form of the full time-evolution operator, e lL ^ t ~ t °\ 
which operates from the left in the second term on the r.h.s. of this equation, must be 
e %L(t-s)^ w \y exe s j s an integral variable. |||| Therefore, this equation is called the TCL 
equation. 

When we ignore the time dependence of the projection operator, Eq. is the same 
as Eq. (2-31) in Ref. ||. It should be noted that not Q(t) but Q(t ) operates from the left 
in the third term on the r.h.s. of the equation. Hence, this term represents the effect of 
the initial correlation of the initial density matrix. This point becomes more clear in the 
next section. In the usual projection operator method, such a term is interpreted as the 
fluctuation force. Therefore, the third term can be regarded as the fluctuation force also in 
this time-dependent projection operator method. 

III. COUPLED TRANSPORT EQUATIONS 

Now, we apply the time-dependent projection operator method to a quantum field theo- 
retical model and derive a transport equation. We consider the Hamiltonian with two boson 
fields 

H = H a + H w + Hj, (30) 

where 

H a = J rf 3 xi{$ 2 (x) + (V0(x)) 2 + ml<p 2 (x)}, (31) 
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H n = J d^ l -{Yi\x) + (Vtt(x)) 2 + m 2 y(x)}, (32) 
Hi = [ d 3 xg7r 2 (x)(p(x). (33) 



Here, and H(x) are the conjugate fields of <f>(x) and ir(x), respectively. We call the 

particle represented by the <f>(x) field the a boson, and that represented by the ir(x) field the 
7r boson. The nonperturbative and interaction Hamiltonians are given by H = H a + H 7T and 
Hi, respectively. Here, we ignore the time dependence of the nonperturbative Hamiltonian 
for simplicity. 

These four fields are expanded as 

0(x,t o ) = E -F=^M*o)e ikx + 4(to)e- ikx ), (34) 



$(x,t ) = -i£ ^(a k (t )e lkx - a| { (to)e- jkx ), (35) 



k 



vr(x,t ) = X: ^^(6 k (to)e ikx + & k (*o)e~ ikx ), (36) 



2VujI 



H(x,t ) = ^E\/^( fo k(to)e lkx - 6U^o)e- lkx ), (37) 

where a; k = ^k 2 + m 2 and = ^/k 2 + m 2 . Here, V and t are the volume of the total 
system and the initial time at which we prepare the initial state. We take the limit V — > oo 
at the end of the calculation. The creation and annihilation operators of the a and tt particles 
are subject to the following commutation relations: 

M*o), 4'(*o)] = [6k(to), bl(t )] = 5 k 3 k ,. (38) 

Here, [ ] represents the commutator. All other commutators vanish. 

For simplicity, we consider the case in which the initial density matrix is given by the 
direct product of the a boson and the n boson density matrices: 

PO = POa ® POtt- (39) 

This means that there is no initial correlation between the a boson and the tc boson. 

The system we are considering consists of two degrees of freedom, that of the a boson and 
that of the n boson. In the usual time-independent projection operator method, it is often 
assumed that one degree of freedom plays the role of a heat bath in thermal equilibrium 
for another degree of freedom. |TJ| In the present application, we treat the two degrees 
of freedom on an equal footing and derive coupled transport equations. [|T0|-|r2"|| First, we 
calculate the transport equation of the a boson. In this case, the a boson is regarded as the 
system and the 7r boson as the environment that is to be coarse-grained, respectively. To 
integrate out the environment degree of freedom, we define the time-dependent projection 
operator as 

P(t)0 = Tr E [p E (t)0}, (40) 
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where 

PE (t)=e-^ H -/Z n (t), (41) 
Z ff (t) = Tr[e-^ H *}. (42) 

Here, Pe{1) is the local equilibrium density matrix. This projection operator satisfies the 
conditions ( |T3"D and ([TU]) . This choice of p_g(t) comes from our implicit assumption that 
the time evolution of the tt boson is well approximated by the local equilibrium density 
matrix with the time-dependent temperature /3~ l (t). It is possible to calculate the transport 
equation without assuming local equilibrium in the time evolution of the system. However, 
in such a case, it is necessary to calculate a large number of correlation functions. To avoid 
this difficulty, we assume local equilibrium in this paper. Now, we employ the following 
initial condition: 

PE(to) = POvr- (43) 

This condition is needed to eliminate the contribution from the third term on the r.h.s. of 
Eq. This becomes clear in the next paragraph. The time dependence of the temperature 
is determined later. 

Substituting O(t ) = a\,(t )ak(to) into Eq. (|i~2"D and using the definition (|40|) , we obtain 



ig 



(pk(Mo))tMt)-4(t)) 

2Vu£ 

+ JTT^ r^{([p-k( S ,to),Pk(t,to)] + )^ (s - t) + ([ Pk ( S ,t ),P-k(t,to)] + )*e^ (s - 4) 
-2<p_ k (s, t ))t(p*(t, to))^-* 5 - 2(p_ k ( S , t )) t (p k (t, t )) t e-^ s -V} 



9 2 « 



^{([P-k(s,to),Pk(t,to)]) i k(t),a_ k (t)e-^^ + a| c (t)e^( s -')] 



Wu£ Jto 

~(W(s, t Q ), p_ k (t, t„)]>t[4(*), a k (t)e"^^ + aUt)e^'-^] + } 

+flu, (44) 

where 

(p k (Mo))t = J rf 3 xe' kx Tr[p £ (t)vr 2 (x,t;t )], (45) 
([p kl Mo),p k2 (Mo)]h = J d 3 x a d 3 x 2 e iklXl e ik3Xa 'ft^(t)[7r 2 (x a ,s;to),7r 2 (x aj t;to)]]. 

(46) 

([p Ul (Mo),Pk 2 (Mo)]+)t = J d 3 x 1 d 3 x 2 e iklxl e ikaxa, n:[p £; (*)[7r 2 (x 1 ,s;to),7r 2 (x 2 ,i;to)] + ]. 

(47) 



Here, [ ] + is the anti-commutator. The last term, flu, on the r.h.s. of Eq. (44]) comes 
from the third term on the r.h.s. of Eq. (£2^). When we calculate the expectation value by 
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using the initial density matrix, this term vanishes because of the definition of flu and the 
condition (|43f) . When we consider the initial density matrix with an initial correlation, we 
cannot employ the condition ([39|), and the third term has a finite value. In short, the third 
term represents the effect of the initial correlation. We can thus regard the third term flu 
as the fluctuation force term as in the usual projection operator method. In this calculation, 
we do not consider the initial correlation, and therefore, we ignore this term. 

Before continuing with the calculation, we discuss the structure of this equation. This 
equation has strange terms that are proportional to a k (£), a^t), aj c (t)a^ k (t) and cik(£)a_k(£). 
Such terms are not seen in the usual Boltzmann equation, which is constituted of terms pro- 
portional to a k (t)dk(i). The terms proportional to a k (t) and cik(t) survive when we consider 
the case that there is a condensate of the a boson. Here, we do not consider the case of the 
condensate, and therefore we ignore such terms. The existence of the terms proportional to 
a k (£)al k (£) and a^{t)a^{t) was pointed out in Ref. fpL5|| . In that paper it is noted that the 



existence of such terms is homologous to a parametric amplifier, and ak(^) a -k(^) 5 a k(£)a-k(t) 
and a k (t)ctk(i) form the SU(1,1) symmetry group. In this calculation, we drop such terms 
for simplicity. 

To calculate the correlation functions ([45]) , (|46|) and ((47D , it is convenient to use a 
technique of thermo- field dynamics (TFD). (The detailed calculation is given in Appendix 
0.) Simply quoting the result, we obtain 

Tr[p E (t)p k (t, t )] = / rf 3 xe ikx 7r(x, t) ^(x, t) (48) 

J I 1 B(t) 

TV[p £ (t)[p kl (s,to),Pk 2 (t,to)]] = / d 3 x 1 d 3 x 2 e iklX1 ^ 



/ 3 Xl d 3 x 2 e iklXl e ik2X2 
x2{(tt(xi, s) tt(x 2 , t)f - (tt(x 2 , t) tt(xi, s) ) 2 }, 

I lB(t) 1 'B(t) 



(49) 



where 



7r(x ll t 1 )7r(xa,t 2 )=£^((l + ^^ 
1 i B (t) k ZVLO^ 



(50) 



Substituting the above results into Eq. (|44j) and taking the expectation value with respect 

ix, we obtain tli 

r <i 1 

ds / cTa 



to the initial density matrix, we obtain the transport equation of the a distribution function, 

d a 9 2 

x [{(1 + nl){l + < +k )(l + <) - nX+k^k} cos« + <+k + ~ t) 

+{« +k (l + <) " (1 + + < + k)<} cosK + ^q+k - Wk)(s - *) 
+{(1 + n>; +k (l + O ~ Ki 1 + <+kK) cos« - < +k + - t) 



+{n*(l + n; +k )(l + K) - (1 + n>£ +k <} cos(u£ - <+k - - «)], 



(51) 



where n k (0 = Tr[pa k (t )cik(*o)]- Here we have omitted the time dependence of n k and n n 



k 



for simplicity. Each term in the braces on the r.h.s. of this equation is composed of two 
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contributions: the gain term and the loss term. The first term describes the creation of two 
tt bosons and one a boson minus their annihilation. The second term describes the creation 
of one a boson and the annihilation of two tt bosons minus the annihilation of one a boson 
and the creation of two tt bosons. The third and fourth terms describe the creation of one a 
boson and one tt boson and the annihilation of one tt boson minus the annihilation of one a 
boson and one tt boson and the creation of one tt boson. This is a reasonable result, which 
is expected from the usual Boltzmann equation. P~7|JT5[| When we take the limit to — > — oo, 



Dirac delta functions that preserve energy conservation are obtained. 

To solve the above transport equation, it is necessary to determine the time dependence 
of the tt distribution function. Our next task is to obtain the transport equation for the tt 
boson. Our strategy is to treat the a and tt degrees of freedom on an equal footing. We 
thus reverse the roles of the a and tt degrees of freedom; that is, we regard the tt boson as 
the system and the a boson as the environment and calculate the transport equation of the 
tt boson using the same procedure as in the calculation of Eq. fl5*T|). In this way, we obtain 
the coupled transport equations of the a and tt distribution functions. 

The projection operator is given by 

P(t)0 = Tr E [p E (t)0}, (52) 

where 

p E (t) = e-M> B '/Z+(t), (53) 

Z <l> {t)=Tr[e- f3 ® H <'}. (54) 

As in the case of the a boson, the a distribution function is included in the tt transport 
equation. We assume local equilibrium in the time evolution again, i.e., 

Tr[p E (t)al(t )a k (t )} =nftt), (55) 



where n^(t) is given by the solution of Eq. fl5i~D 

When we use the definition (|53"|) and substitute O(t ) = tf(to)b(to) into Eq. ([12]) , we can 



obtain the transport equation of the tt distribution function. In this case, the transport 
equation has fourth order correlation functions of the tt boson. To obtain a closed form for 
the coupled transport equations, we approximate such terms as 

Tr[p bt(t)bl(t)b m (t)b n (t)} = ^(tK(t)Og +nmnUW£l&- 



(56) 



Finally, the transport equation for the tt distribution function is 
d -nl= tds fdH- ^ 



dt k Jt J (2ir) 3 u^uf +ll 

x [{(1 + < + i)(l + Oil + <) - cosK +k + «T + 0(s - t) 

+K + i(l + K)(l + <) - (1 + <+,)«} cosK +k - wf - c£)(* - t) 
+{(1 + < +1 )(l + <)< " + <)<} cosK +k - + wD(a " t) 

+K+i(l + K)n\ - (1 + + nfK} cosK +k + wT - - t)]. 



(57) 
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Each term in the braces on the r.h.s. of Eq. (|57 ) can be interpreted as gain minus loss 
processes, as in the case of the a transport equation. The time-evolution of the a and tt 
distribution functions are obtained by solving this coupled set of equations. 



IV. SUMMARY AND CONCLUSIONS 

We have derived a systematic perturbative expansion formula for a Langevin-type equa- 
tion without time-convolution integral terms. In this formalism, the irrelevant subsystem 
can have some time dependence, and therefore, we can treat a more general nonequilibrium 
process in which the mass, the temperature, and so on, are time dependent. When we 
ignore the time dependence of the projection operator, we can reproduce the result of the 
usual projection operator method. |J Furthermore, the third term on the r.h.s. of Eq.(p9|) 
represents the effect of the initial correlation. We can thus interpret it as the fluctuation 
force even in the time-dependent projection operator method. 

We applied this formalism to a quantum field theoretical model that consists of a and 7r 
bosons and thereby obtained coupled transport equations for the two bosons. Deriving these 
coupled transport equations, the correlation functions were calculated with respect to the 
local equilibrium density matrix. In other words, we have assumed that the time evolution 
of the system can be well approximated by the local equilibrium distribution function. The 
derived equations have terms that are not seen in the usual Boltzmann equation. It is worth 
studying the effect of such terms on the transport equation, but we have dropped them in 
this paper. Each transport equation has a form that can be interpreted as gain minus loss 
processes, which is the structure usually seen in the Boltzmann equation. This gives reason 
to believe that the time-dependent projection operator method developed in this paper is a 
valid formalism to describe nonequilibrium processes. 

In this paper, we have ignored the time dependence of the nonperturbative Hamiltonian. 
In general nonequilibrium processes, it is possible to consider the situation in which the mass 
of a particle changes with time. In this case, the time dependence of the nonperturbative 
Hamiltonian is caused by the time- dependent mass. We will report on an investigation of 
this effect in a future publication. 

Finally, we mention the fluctuation-dissipation theorem. This is a famous theorem in sta- 
tistical mechanics that expresses the relation between macroscopic transport coefficients and 
microscopic fluctuations. It is well-known that there are two expressions for the fluctuation- 
dissipation theorem. For the Langevin equation, the fluctuation-dissipation theorem of the 
second kind is important. This theorem was first proven by Mori using a time-independent 
projection operator technique. Furukawa developed the time-dependent Mori projection 
operator, according to which the generalized version of the fluctuation-dissipation theorem 



of the second kind was obtained. [14]] The time-dependent Mori projection operator, which 



satisfies P(t)P(t') = P(t), is not included in our formalism, due to condition (|T3]). Our 
formalism includes the time-independent projection operator method and reproduces the 
Tokuyama-Mori equation H when we use the Mori projection operator. Therefore, the 
fluctuation force in our formalism agrees with that in the Tokuyama-Mori formalism. We 
thus conclude that our fluctuation force is a natural extension of theirs, which satisfies the 
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fluctuation-dissipation theorem. However, the fluctuation-dissipation theorem in our formal- 
ism has still not been established, and demonstrating it is one important future problem. 
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APPENDIX A: DEFINITION OF THE OPERATORS C AND V 

The definitions of the operators C(t,to) and V(t,t ) are 

C(t,t ) = Uo(t,t )^ tLit ' to \ (Al) 

. dsLQis) ^ 

v [h t )=e^> (U$)-\t,t ), (A2) 



where 



i f* dsLo(s,to) 

[/ (Mo)=e> , (A3) 

„ i f* dsL (s,t )Q(s) 

£#(f,<o)= e > ' (A4) 

These operators must satisfy the following differential equations: 

jC{tM = U (t,t )(iL (t,t )-iL)e- iL ^ 

= -iLj(Mo)C(Mo), (A5) 

A i f* dsLQ(s) „ i 

-P(t,t ) = e> (<L -*L (Mo))QW(t/ V(Mo) 



V(t,t )iL?(t,t ), (A6) 



where 



Li{t,t ) = ^(Mo^KMoK^Mo), (A7) 
L?(Mo) = £/ °(t, t )M*> *o)<5(^o Q ) _1 (^ to)- (A8) 

The solutions of the above differential equations are 

C(Mo) = 1 + EH)" / ^1/ ^•••/ dt n L I (t 1 ,t )L I (t 2 ,t )---L I (t n ,t ), 

(A9) 

£>(Mo) = l + Y J i n f dt! r dt 2 --- I*'' dt n Lf(t 1 ,t )Lf(t 2 ,t )---Lf(t n ,t ). 

n=l 

(A10) 
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APPENDIX B: TRANSFORMATION OF THE OPERATOR £(T,T ) 
The operator £(i, i ) can be rewritten as 

i J^drLQ(r) 



£(Mo)= / dse- iL{t - s) {Q(s) + P(s)iLQ(s)}e 
J to 

J t as Jt as 

= -P(t) + e-^-^P(t )S dTLQ(T) + f dse-^- s hLP(s)eS d " 

Jtn 



drLQ(r) 

Vltj-e — ""WojeV 



to 



(is 



z / dsLo(s,to)Q(s) 

= Q(t) - U Q l {t,t )C{t,t )Q{t )V(t,t )ed ta . (Bl) 

Using mathematical induction, we confirm the following relation: 

P(t)£(Mo)(QW£(Mo)) n 



[{-l)^P{t){Q(t)(C(t) - l)TQ{t) + 


r- l p(t){{c{t) - i)Q(0} n+1 


n-l 
1=0 


-iXq^MMo))"- 1 -' 


n-l 

+^(*) £(-!)'{(<?(*) - i)g(t)} m (^(t) - 

2=0 


llWfimt.to))- 1 - 1 


-p(*)E(-i)'{g(*)(c(t) - i)}'Q(t)(x>(t) 


-l)(Q(t)£(Mo)) n ^ 


-^) E(-l) ; {(C(t) - l)Q(t)}< +1 (X>(*) - 
z=o 


l)(Q(t)E(t,t )) n - z , 



(B2) 

z=o 
where 

Q(t) = U^(t,t )Q(t )U^(t,t ), (B3) 
= IV^MoWMoWMo), (B4) 
= (V?)-%h)V(t,h)Vt?(t,to). (B5) 

Here, n is an integer and n > 1. In this derivation, we have used the relation 

Q(*)Q(t) = Q(t), (B6) 

which can be proved from the condition (PHI). The second and third terms in 
P(t)S(t,t )(g(t)S(t,t )) n and the fourth and fifth terms in P(t)E(t, t )(Q(t)E(t, to))"" 1 
cancel. The fourth and fifth terms in P(i)£(£, to)(Q(t)H(t, to)) n and the second and third 
terms in P(£)E(t, to)(Q(t)H(t, to)) n+1 also cancel. Therefore, only the first term survives. 
As a result, all the terms including T>(t) disappear. 
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Note that the relation E(t, to) = T,(t, to)Q(t) can be derived from Eq. (p!6[). We thus find 
= P(t)E(t,t ): ' 



-Q(i)£(Mo) 

oo 

p(*)£(Mo) £(£(*)£(*, *o)) n 



n=0 

oo 



-P(t) £[{-Q(*)(C(*) - VVQW - {-(C(t,t ) - i)Q(t)} n 

n=0 

oo 

-P(t)U \t,t ) £[{-Q(t)(C(t,t ) - l)} n Q(t) 



n=0 

i n+li 



-{-(C(i,«„)-l)Q(i)}" +1 ]Co(«,«o) 
-F M[ > 1 -( t , il )«t) 1 + (c(Mj ' ) _ lW(i) g,(M.) 



-P( i)t / -'(Mo)(C(M„) - ')QW 1 + (c(f , f , 1 ) _ 1)0W W.t 

1 



P{t)UM t t M t, t )Q(t) 1 + _ Dbft, to). 



(B7) 



APPENDIX C: CALCULATION OF CORRELATION FUNCTIONS BASED ON 

NONEQUILIBRIUM TFD 



To calculate the correlation functions (|45|) , (|46f) and (|47D , it is convenient to use a 
technique of thermo-field dynamics (TFD). In TFD, the statistical average is expressed as 
a kind of a vacuum expectation value. As an example, we consider the statistical average of 
the system with the free Hamiltonian of a boson system, H = J2 cu k d k d k . First, we introduce 
the transformation through which the new pairs of creation and annihilation operators D k 
and D k , and and _D k are defined: 

4 = cosh0 k L>t + sinh0t.Dk, (CI) 

d k = cosh6 k D k + sinhflkDj, (C2) 

where 

sinh 2 ^ = n k , cosh 2 ^ = 1 + n k . (C3) 

Here, n k is a Bose distribution function. The above operators satisfy the commutation 
relation 

[D k ,Dl] = [D k ,Dl} = 6gl,. (C4) 
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All other commutators vanish. Now, we can define the thermal vacuum \9d) as 

D±\6 D ) = D k \9 D ) = 0. (C5) 

We thus express the statistical average in terms of the vacuum expectation value: 

Tv[p th O] = (0 D \O\0 D ), (C6) 

where 

pth = e-W/Z, (C7) 
Z = Tr e~ m . (C8) 

Next, we return to the calculation of the correlation functions with the local equilibrium 
density matrix. In the case of a time-dependent temperature, we introduce a time-dependent 
transformation to define the new pairs of creation and annihilation operators B^t) and 
Bl(t), and B k (t) and £?£(t): 

&£(*„) = cosh6Z(t)Bl(t) + S inh0l(t)B k (t), (C9) 

& k (t ) = cosh9l(t)B k (t) + sinh0£(t)5 k (t), (CIO) 

where 

sinh 2 0J(t) = nl(t), (Cll) 

cosh 2 0£(t) = l + nl(t). (C12) 

Here, n£(t) is the time-dependent distribution function of the 7r boson: 

nl(t) = T^(*)&L(*o)&k(*o)] 

' (C13) 



e 0(*K - 1 ' 

The creation and annihilation operators B k (t), B k (t), B k {t) and B k (t) satisfy 

[B k (t),Bl(t)} = [B k (t), Bl(t)} = Sl%. (C14) 
All other commutators vanish. Then, we can define the time- dependent vacuum |#b(£)) as 

B k (t)\e B (t)} = B k (t)\e B (t)} = o. (ci5) 

We thus express the statistical average with the local equilibrium density matrix in terms 
of the time-dependent vacuum expectation value: 

Tr[p E (t)0] = (9 B (t)\0\6 B (t)), (C16) 

where pE{t) is defined in Eq. (|4lD . Therefore, we can simplify the calculation of the corre- 
lation functions with the help of Wick's theorem. Simply quoting the result, we obtain 
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Ti[p E (t)p k (t,t )] = (M*)|pk(Mo)IM*)) 

= J rf 3 xe* kx 7r(x,t)^ (x,t) 
Tr[p E (t)[p kl (s, t Q ), p k2 (t, t )]] = (0 B (t)| [/>i(s, *o), Pk(*, to)] IM*)> 

= /" d 3 Xl d 3 x 2 e iklXl e ik2X2 

x2{(7r(x x , s) 7r(x 2 , t)) 2 - (tt(x 2 , t) n(x u s) ) 2 }, 



(C17) 



J B(t) 



J B(t) 



(C18) 



where 

7r(xi,ti)7r(x2,t 2 ) 



'B(t) 



V J_ 



:[cosh^ l (t) J B kl (t)e 



-iug (ti— to) 



+ sinhr kl (t)S_ kl (t)e tok i 



mj£ (tl— to) 



sinh ^ k2 (t)Bl k2 (t)e- l <2 + cosh k2 (t)B^ (t)e^ ( * 2 " t()) ] 
= ^ -^(cosh 2 9l(t)e^<^ t2) + sinh 2 ^(t)e^ (tl -* 2) )e ik(xi - X2) . 

Readers familiar with nonequilibrium thermo-field dynamics may notice the difference be- 
tween the technique that we have outlined here and that used in Ref. [17fl. However, these 
two techniques give the same result with respect to the trace calculation. 



(C19) 
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